Detection of the phase shift from a single quantized superconducting vortex 



Taras Golod, Andreas Rydh & Vladimir M. Krasnov* 
Department of Physics, Stockholm University, AlbaNova University Center, SE - 106 91 Stockholm, Sweden 

An Abrikosov vortex in a superconductor carries a flux quantum, €>o = hc/2e, localized at its 
center, but induces a global 2tt phase rotation in the superconducting condensate. This long-range 
gauge field[l] outside the area pierced by a magnetic field is due to the Aharonov-Bohm effect [2], 
which is a non-classical phenomenon that illustrates the significance of potentials rather than forces 
in quantum mechanics. In the London gauge, the phase of the condensate is given by the polar angle 
around the vortex. Here we raise the question whether this phase shift could be detected by means 
of Cooper pair interferometry using Josephson junctions as phase-sensitive detectors. We introduce 
a single Abrikosov vortex into a superconducting lead with a detector junction made at the edge of 
the lead. We observe that the vortex induces a Josephson phase shift equal to the polar angle of 
the vortex within the junction length. When the vortex is close to the junction it induces a 7r-step 
in the Josephson phase difference, leading to a controllable and reversible switching of the junction 
into the - 7r state[3-12]. This in turn results in an unusual <&o/2 quantization of the flux in the 
junction. The vortex may hence act as a tunable "phase battery" for quantum electronics. 



A sketch of our experiment is shown in Fig. la. A sin- 
gle Abrikosov vortex is placed in a superconducting lead. 
A detector Josephson junction is made at the edge of the 
lead. The supercurrent I c through the junction is a result 
of interference of Cooper pair wave functions, which leads 
to a Fraunhofer modulation of I c as a function of mag- 
netic field. The Josephson phase shift, induced by the 
vortex, is detected from a comparison of I C {H) patterns 
with and without the vortex. 

In the London gauge, variation of the phase of the 
superconducting condensate around the vortex is given 
by the polar angle 0„, which, at the junction interface is 
equal to: 

@ v (x) — arctan ^- — —^j + const., (1) 

where x v and z v are the vortex coordinates and x is the 
position along the junction length. Profiles of @ v (x) for 
different distances from the vortex to the junction arc 
shown in Fig. lb. Even in quantum mechanics gauge 
fields have limited physical significance. Only closed 
path integrals of gauge fields arc measurable [1]. For 
Cooper pairs such integrals around the vortex are equal 
to 2n, which is indistinguishable from in the absence of 
the vortex. Open path integrals are not gauge-invariant 
and should not be measurable. Therefore, the question 
whether a distant Abrikosov vortex gives rise to a Joseph- 
son junction phase shift is non-trivial. 

The main challenge for the present experiment is to 
avoid vortex intrusion into the junction area, which might 
induce a parasitic phase shift in the junction. Conven- 
tional Josephson junctions are formed by a barrier sand- 
wiched between thin superconducting films. Abrikosov 
vortices in such "overlap" junctions tend to minimize 
their energy by orienting themselves perpendicular to the 
electrodes, thus introducing a segment of a Josephson 
vortex (fluxon) in the junction. This is a well-known 
reason for distortion of I C (H) patterns in conventional 



overlap junctions subjected to out-of-plane fields[13, 14]. 
To avoid fluxon formation we employ two types of spe- 
cially designed detector junctions: 

1. Planar Nb/CuNi/Nb junctions, see Fig. Ic. Such 
junctions are ideal for the planned experiment: due 
to their two-dimensional geometry, the Abrikosov 
vortex, which is oriented perpendicular to the Nb 
film, cannot cross the junction line. 

2. Mesoscopic Nb/PtNi/Nb junctions, see Fig. Id. 
Mesoscopic sizes of the junctions confine the vor- 
tex in the middle of the electrode parallel to the 
junction plane[15] and allow the detection of I c in 
very strong magnetic fields (up to 20 kOe), which 
further helps to align the vortex. 

In Fig. 2a we show I c for a planar Nb/CuNi/Nb junc- 
tion as a function of the magnetic flux $ through the 
junction (field is oriented perpendicular to electrodes, 
along the y-axis in Fig. la). Measurements were done by 
first sweeping the field from to 40 Oe, then to —40 Oe 
and finally back to 0. The J r c ($) patterns are almost iden- 
tical for all three sweeps, except for an offset A$, which 
changes stepwise with the field, as shown in Fig. 2b. The 
apparent quantization of A$ implies that each step is 
caused by entrance or removal of an Abrikosov vortex in 
the electrodes. Remarkably, the offset A$ is quantized 
in half flux quanta. As a result, the I c ($) modulation 
gets out-of-phase, i.e., positions of minima and maxima 
are interchanged, with each step in A<F Furthermore, 
the offset occurs in the direction of applied field, which 
means that the trapped field in the junction is opposite 
to the direction of the field. 

To clarify the origin of the unusual $o/2 quantiza- 
tion, the junctions were modified by a Focused Ion Beam 
(FIB) in two steps, as shown in Fig. Ic. First, a vor- 
tex trap was made in order to control the position of the 
vortex. The trap is a small hole ~ 30 nm in the cen- 
ter of one of the electrodes near the junction. Second, 



FIG. 1. Experimental configuration and sample geometry, a, 

Geometry of the experiment: the phase shift from the single 
Abrikosov vortex is detected by the junction at the edge of the 
lead. The polar angle of the vortex within the junction length 
AQ V is marked by the dashed lines, b, The polar angle of the 
vortex Q v (x) along the junction length for different distances 
z v from the vortex to the junction and x v = L/2. c, Top view 
of a planar Nb/CuNi/Nb junction #1, with a vortex hole 
and a stray field drain, d, SEM image of a nano-sculptured 
Nb/PtNi/Nb junction. SEM images in c and d are shown in 
the same perspective as the sketch in a. The magnetic field in 
our experiment is applied along the y-axis (into the paper). 



a stray field drain was added by removing a substan- 
tial part of the electrode in the vicinity of the vortex 
hole. The drain should substantially decrease the mag- 
netostatic stray fields from the vortex at the junction. 
Figure 3a shows the I C (H) patterns after these modifi- 
cations in the absence of Abrikosov vortices. It is seen 
that the shape of the I C {H) pattern remained intact af- 
ter introducing the field drain, implying that the junction 
uniformity was not affected by the drain. 

A vortex in the hole can be controllably introduced or 
removed by applying an appropriate field and by using 
Lorentz force from the transport current (see supplemen- 
tary information). Figure 3b shows a new type of I C (H) 
patterns, which appears after trapping the anti-vortex in 
the hole. It has three new characteristic features: 



1. The central maximum at H 
minimum. 



is replaced by a 



2. The I C (H) modulation is out-of-phase with that for 
the vortex-free pattern, i.e., A$ ~ 3>o/2. 

3. The periodicity of the I C {H) modulation doubles at 
the left side of the pattern, leading to a clear left- 
right asymmetry [16]. When a vortex is trapped 
in the hole (instead of an anti- vortex) , the I C (H) 
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FIG. 2. Quantized flux offset in Fraunhofer patterns I c ($) 
of a planar junction without vortex hole or stray field drain. 

a, measured patterns for consecutive field sweeps from 

to 40 Oe (red), from 40 to -40 Oe (blue) and from -40 
to Oe (green). The appearance of hysteresis (flux offset) 
upon sweeping the field is clearly seen. Note that the / c ($) 
modulation for the first sweep (red curve) is out-of-phase with 
that on the way back (blue curve), which indicates that the 
offset is half-integer of $o- b, Measured flux offset vs. H(&) 
for the same field loop (hatched area corresponds to the range 
shown in panel a). Symbols represent minima or maxima in 
I c ($). The unusual 4>o/2 quantization of A$ is clearly seen. 
Each step corresponds to a sequential entrance or exit of one 
Abrikosov vortex. The expected vortex configurations are 
indicated by adjacent sketches. Note that the offset of I c {&) 
patterns occurs in the direction of the applied field, which 
implies that the trapped flux in the junction is opposite to 
the applied field. 



pattern becomes mirror reflected with respect to 
the H = axis (not shown). 

These are the well known fingerprints of - tt junctions, 
with a step-like 7r-shift in the Joscphson phase difference 
within the junction[3-6]. Properties of tt -junctions with 
negative Josephson coupling and - tt junctions have 
attracted significant attention in recent years, both due 
to the interesting physics involved, and the potential for 
new applications. So far, three types of - tt junc- 
tions were realized based on: i) the d-wave symmetry 
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FIG. 3. I C (H) modulation for a planar junction with the 
vortex hole and the stray field drain. Sketches demonstrate 
junction sample geometries: blue curves - with the vortex hole 
only, red curves - with the hole and the stray-field drain. The 
current is normalized on the maximum critical current in the 
absence of vortices, I c o- a, Without vortices, b, With an anti- 
vortex in the hole. To trap the anti- vortex, a sufficiently large 
negative field was applied before the measurements. Clear 
signatures of the - 7r state are seen: i) the central maximum 
is replaced by a minimum, ii) modulation of I C (H) gets out- 
of-phase with that without the vortex and iii) doubling of 
the periodicity occurs at one side of the pattern. Note that 
introduction of the stray field drain affects neither the I C {H) 
pattern, nor A$. 
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FIG. 4. Influence of the junction geometry on the flux offset 
in mesoscopic Nb/PtNi/Nb junctions, a, b, I C (H) patterns 
for the same Nb/PtNi/Nb junction 1140 x 230 nm 2 for two 
in-plane field orientations. First the field is swept up starting 
from the vortex free state until first Abrikosov vortices en- 
ter the electrodes (blue symbols) and then back to zero (red 
symbols). Entrance and exit of the Abrikosov vortex leads 
to a sudden appearance of the offset A$ in the Fraunhofer 
pattern. The inset in panel b shows the I C (H) pattern in a 
wider field range, c, Flux offset A$/$o induced by the vortex 
as a function of the polar angle A6„ of the vortex within the 
junction. Similar symbols correspond to the same junction 
with different in-plane field orientations. Red squares corre- 
spond to the junction shown in panels a and b. The dashed 
line indicates that the Josephson phase shift induced by the 
vortex (right axis) is simply equal to the polar angle of the 
vortex. 



of the order parameter in high T c superconductors [7- 
9], ii) the oscillatory nature of the proximity induced 
order parameter in hybrid superconductor/ferromagnet 
junctions[5, 6, 10, 11], and iii) the phase shift by current 
injection into the junction[4, 12]. Here we demonstrate 
that a conventional 0-junction can be switched into the 
- 7r state by a single Abrikosov vortex, parallel to the 
junction plane. 

Importantly, introduction of the stray field drain does 
not reduce the offset A<I>, see Fig. 3b. This clearly shows 
that the effective trapped flux in the junction is not 
driven by the simple magnetostatic spreading of the vor- 
tex stray field. To get more insight into the influence 
of junction geometry on A<3>, mesoscopic Nb/PtNi/Nb 
junctions were used, where the vortex is geometrically 
confined in the middle of the electrodes [15], z v ~ d/2, 
x v ~ L/2. All junctions had the same electrode thick- 
nesses d, but different junction lengths L, thus allowing 



variation of the ratio z v /L ~ d/2L. 

In Figs. 4a,b we show I C (H) for a Nb/PtNi/Nb junc- 
tion 230 x 1140 nm 2 with an in- plane field parallel to 
different facets of the junction. The A$ induced by the 
vortex is clearly seen in both cases. For the long junction 
case, L = 1140 nm, A$ ~ $0/2, which results in out-of- 
phase modulation of the patterns. However, for the short 
junction case, L — 230 nm, A$ ~ 0.16<I>o is considerably 
smaller. 

Figure 4c summarizes the measured A$ for all stud- 
ied Nb/PtNi/Nb junctions. Obviously, A<t>, introduced 
by the vortex increases considerably with the junction 
length L. The observed offset A$ in I C (H) corresponds 
to the net variation of the Josephson phase difference 
along the junction length 

A^(L-O) =2ttA$/$ . (2) 
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Thus the Abrikosov vortex does induce a measurable 
Josephson phase difference in the junction. The dashed 
line in Fig. 4c indicates that the Josephson phase shift 
induced by the vortex is agreeing well with the polar an- 
gle of the vortex within the junction: Aip v (L — 0) ~ A6„ 
(see Fig. la and equation (1)). The A6„ was calculated 
by assuming that the vortex is placed in the middle of 
an electrode z v ~ d/2, where d ~ 300 nm is the average 
thickness of Nb electrodes. 

To understand the origin of the observed phenomenon 
we first exclude several unsustainable scenarios: 

1. The vortex docs not introduce a segment of a 
Josephson fluxon in the junction, which might dis- 
tort I C (H)[13, 14]. In our structures, the Abrikosov 
vortices were oriented strictly parallel to the junc- 
tion interface and never cross the junction area. In 
particular, such crossing is impossible for the planar 
junctions, because of the two-dimensional junction 
geometry. 

2. A direct field of the vortex stretching into the junc- 
tion or magnetism in the junction cannot explain 
the experiments, because the induced field in the 
junction is opposite to the applied field. Indeed, 
from Fig. 2a it is seen that the central maximum 
in J c ($) is shifted to a positive field after apply- 
ing +40 Oe and to a negative field after applying 
—40 Oe. Since the central maximum corresponds 
to zero total flux in the junction, the additional 
field within the junction is always opposite to the 
applied field. 

3. The magnitude of the signal can hardly be ex- 
plained by a finite vortex current at the junction 
interface. Although it produces a phase shift of the 
proper sign, its magnitude should decay strongly 
with the distance from the vortex to the junction[4] 
and should for no reason produce a quantized <&o/2 
flux offset in the junction. Similarly, it is not pos- 
sible to explain the characteristic equality between 
the Josephson phase shift and the polar angle of 
the vortex, shown in Fig. 4c. 

4. Magnetostatic stray field from the vortex would 
also give a phase shift with correct sign. The exper- 
iment with the stray field drain, however, demon- 
strates that the Josephson phase shift is unaffected 
by variation of magnetostatic conditions. Further- 
more, there is no reason for the stray field to be 
quantized as &o/2. When the vortex is placed very 
close to the junction, clear signatures of the - tt 
junction are seen, see Fig. 3b. This implies that the 
induced Josephson phase shift has a form of a sharp 
tt step, which is again difficult to explain in terms 
of simple magnetostatics because it would require 
field focusing in one point. Additional discussion 
can be found in the supplementary. 



Our data show an unambiguous correspondence be- 
tween the Josephson phase shift and the polar angle 0„, 
which represents the variation of the phase of the su- 
perconducting condensate around the vortex within the 
London gauge, equation (1), as if the phase of the conden- 
sate is rigidly coupled to rotation of the current in the 
vortex. Numerical simulation presented in the supple- 
mentary demonstrate that equation (1) provides a good 
overall agreement with all our observations. It naturally 
explains the unusual half-integer flux quantization in the 
junction. The associated tt- Josephson phase shift in this 
case is simply equal to the change in the polar angle upon 
going from the left to the right side of the Abrikosov vor- 
tex. When the vortex is close to the junction, the po- 
lar angle changes stepwise, as shown in Fig. lb, and the 
junction switches into the - tt state. Yet, note that 
the remarkable success of the London gauge description 
of the phase shift around the vortex is surprising because 
phase shifts between any two points (such as the left and 
right edges of the junction) are not gauge invariant and, 
therefore, should not be measurable. We assume that 
the presence of the junction as such plays a crucial role 
on the way from the unmeasurable phase shift of the su- 
perconducting condensate to the measurable Josephson 
phase difference. Although the seeming rigidity of the 
gauge field around the Abrikosov vortex remains to be 
clarified, we demonstrate that it can be employed as a 
tunable and reversible phase battery for Josephson elec- 
tronics. Depending on the geometrical factor z v /L, such 
a battery can provide either a quantized step-like 7r-shift, 
or an arbitrary phase shift in the range < Aip < tt. 

METHODS 

Planar Nb/CuNi/Nb junctions of the "variable thick- 
ness" type were made by cutting CuNi/Nb double layers 
by a Focused Ion Beam (FIB). The thicknesses of CuNi 
and Nb layers were 50 and 70 nm respectively. The junc- 
tion length and width were ~ 4 /im and 20 nm respec- 
tively. The details of junction fabrication and charac- 
terization can be found in Refs. [18, 19]. Mesoscopic 
Nb/Pti_a;Ni x /Nb junctions of the "overlap" type were 
made by 3D FIB-sculpturing. The thicknesses of bottom 
and top Nb layers were 225 and 350 nm respectively. 
The thickness of the Pti-zNLj layer was varied from 20 
to 24 nm and the junction length was varied from 110 to 
1200 nm. To control I c , the Ni concentration was varied 
from to 67 at.%. In Fig. 4c we show data for junctions 
with pure Pt barrier (circles and down triangle) and with 
54 at.% (squares), 60 at.% (up triangle) and 67 at.% (dia- 
mond and star) of Ni. Details of junction fabrication can 
be found in Ref. [19]. Nb, CuNi and PtNi layers were de- 
posited by magnetron sputtering on oxidized Si wafers. 
The films were first patterned into 5 /im wide bridges 
by photolithography and ion etching (CF4 reactive ion 
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etching for Nb and Ar- milling for CuNi and PtNi). 

Measurements were done in a four probe configuration 
in liquid 4 He or flowing gas cryostats at T ~ 1.8 K. I c 
was measured by ramping the bias current and detecting 
the largest current corresponding to a voltage less than 
a certain threshold level. 
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SUPPLEMENTARY INFORMATION 

Controllable manipulation of the vortex by 
transport current 

Transport current exerts a Lorentz force on a trapped 
vortex. This allows controllable manipulation of the vor- 
tex state[l]. We can remove a vortex by applying a large 
enough current at zero magnetic held, or introduce a vor- 
tex by applying current at finite magnetic field. The main 
panel in Fig. 5 shows I C (H) for the same Nb/PtNi/Nb 
junction as in Fig. 4b, measured at two maximum bias 
currents lb (the current was swept in a saw-tooth like 
manner from to > J c ). The current lb = 0.5 mA 
is too small to cause vortex entrance and the system 
stays vortex-free in the shown field range (solid symbols) . 
However I b = 2.2 mA is sufficiently large to manipulate 
the vortex state and at H ~ 0.7 kOe we can switch be- 
tween the vortex-free and the single Abrikosov vortex 
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FIG. 5. I C (H) modulations of a mesoscopic junction of 
size 230 x 1140 nm 2 at T = 1.8 K for bias currents 

I b = 0.5 mA and 2.2 mA. The magnetic field was applied 
across the 1140 nm junction facet and swept up. Inset. I c (h) 
dependence at H — 0.69 kOe. The bias current shakes the 
system to change the vortex state. By suddenly decreasing 
the current, the vortex state freezes and I c can be measured. 



state solely by the bias current, as indicated by the ar- 
row in Fig. 5. 

Furthermore, we can controllably introduce or remove 
the vortex by changing the maximum bias current. The 
inset in Fig. 5 shows I c as a function of lb at constant 
in-plane magnetic field H — 0.69 kOe. It is seen that by 
changing If, from ~ 2.1 to ~ 2.3 mA the system can be 
switched back and forth between the two vortex states, 
characterized by different I c . Our mesoscopic junctions 
thus act as memory cells (which can be considered as non- 
volatile due to the quantized nature of the vortex). The 
possibility to control the induced flux offset A$ in the 
detector junction solely by the transport current clearly 
demonstrates that the observed phenomenon is caused 
by the Abrikosov vortex. 



Analysis of I C (H) modulation 

Assume that a single Abrikosov vortex is placed in the 
electrode 1 parallel to the junction plane, as shown in 
Fig. la. Supercurrent densities in electrodes 1 and 2 are 
determined by the second London equation: 



A 2 = 



4ttA? i2 



2tt 



^01,-3 



(3) 



Here and A are the corresponding scalar (phase) 
and vector potentials, and Ai^ are London penetration 
depths of the electrodes. The gauge-invariant Josephson 



phase difference ip is obtained by integrating equation (3) 
over an infinitesimal contour of length dx, covering the 
barrier: 
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dx 



2ir 
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(4) 



Here Jf 2 ' are the x-components of the supercurrent den- 
sity in the vicinity of the barrier in electrodes 1 and 2, B 
is the in-plane (j/-axis) magnetic induction in the barrier 
and t is the barrier thickness. 

In the Meissner state, J{ 2 are obtained by solving 
equation (3) in the electrodes, with boundary conditions 
B = H outside the electrodes, z = ±di j2 , where & re 
thicknesses of electrodes. Straightforward calculations 
yield[2]: 



M 



dip 

dx 



2tt 
$ 



BA - HS] 



(5) 



t + X)i=i,2 A, coth f; and S = 
. From comparison of equation (4) and 



Here A = 

£i=i, 2 Aicosechf; 
(5) it follows that current terms in the right-hand side of 
equation (4) represent the effective magnetic flux in the 
electrodes, which depends on the electrode geometry, 
and, therefore, is not quantized. For junctions with bulk 
electrodes di 2 3> A (realized in our planar junctions) 
the H term in equation (5) vanishes and one obtains a 
familiar expression dipm/dx = (2ir/$ )B(t + Ai + A 2 ). 
For e?i j2 <C A, screening by electrodes is weak and 
dip M /dx ~ (27r/$ )B(t + dt/2 + d 2 /2). Integrating 
dip/dx over the junction length x € [0, L], we obtain a 
simple relation between the effective flux and the total 
phase difference shift A<p = p{L) — ip(0) in the junction: 



$ _ Aip 
~ ~ 2tt"' 



(6) 



where <& = d* J Q Bdx is the effective magnetic flux in the 
junction and d* is the effective magnetic thickness of the 
junction, d* = t+X 1 +X 2 for bulk and d* = t+d 1 /2+d 2 /2 
for thin electrodes, respectively [2]. 

Since the London equation (3) is linear, the Josephson 
phase difference tp in the presence of an Abrikosov vor- 
tex in the electrode is simply given by a superposition 
of contributions from the Meissner state, pm(H), equa- 
tion (5), and from the vortex at zero external magnetic 
field, tp v (0): 

p = p M (H) + ip v (0) = ip M {H) + O vl - Q v2 , (7) 

where O^i^ are phases of the superconducting conden- 
sate in each of the electrodes, induced solely by the vor- 
tex. Within the London gauge, 0„i is simply given by 
the polar angle of the vortex, see Fig. la: 



Q v i(x) = arctan 



const. 



(8) 
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respect to the arbitrary constant. The additional con- 
tribution from the anti-vortex to the Josephson phase 
difference was taken to be equal to the polar angle of 
the anti- vortex, equation (8), which are shown in the 
insets. When the anti-vortex is far away from the junc- 
tion, z v ^> L, the I c (&) follows the conventional Fraun- 
hofer pattern with maxima at half- integer "J*/^; as indi- 
cated by the arrow in Fig. 6a. At intermediate distances 
z v ~ L, the anti-vortex introduces a flux offset to I c (^)i 
as indicated by the arrows in Figs. 6b, c. Note that the 
sign of A$ is negative, i.e., it is the same as for the intro- 
duced anti-vortex, consistent with experimental data in 
Figs. 2 - 4 in the main manuscript. At z v <C L, (f v turns 
into a 7r-step and clear fingerprints of the - 7r junction 
appear in the / c ($), which were discussed in connection 
with Fig. 3b in the main manuscript. 



FIG. 6. Simulated ic(3?) patterns with a single anti-vortex 
at different distances from the junction z„ and x v = L/2, 
assuming that the vortex induces the Josephson phase shift 
equal to the polar angle of the vortex. Clear signatures of the 
- 7r junction are seen for z v <C L: The I c {&) modulation 
becomes out-of-phase with the Fraunhofer pattern without 
vortex A$ ~ $o/2, as indicated by vertical arrows in panel c. 
Central minimum appears at $ = in panel d. Doubling of 
the periodicity occurs at one side of the pattern $ < 0. Note 
that the phase gradient (induced field) has the opposite sign 
with respect to the vortex. Insets, Vortex contribution to the 
Josephson phase shift corresponding to the polar angle of the 
vortex. For z v <C L it turns into a sharp 7r-step, which causes 
switching of the junction into the - tt state. 



The phase shift in the second electrode 0„2 is due to 
magnetic flux from the vortex in the electrode- 1, reach- 
ing the electrode-2. Assuming that 9„2 is negligible, the 
Josephson phase difference caused by the vortex, <p V} is 
given the polar angle, equation (8), shown in Fig lb and 
in insets of Fig. 6 for different distances z v from the vor- 
tex to the junction. 

The Josephson critical current in a magnetic field H is 
given by: 



Ic(H) 



IcO 



sin [ip (x)] dx, 



(9) 



where I c q is the maximum critical current at zero field. 
In the Meissner state, v>m(%) is varying linearly, with the 
slope proportional to the applied field, see equation (5): 
Pm{x) — 2ird*B/& + C, where C is an arbitrary con- 
stant. The I C {H) pattern is obtain by maximization 
of equation (9) with respect to C. This leads to the 
Fraunhofer modulation I c ,m{*&) = Ico\ sin(/)|//, where 
/ = 7r$/<I>o, shown in Fig. 6a. 

Figures 6b-d show I c (&) patterns in the presence of 
a single Abrikosov anti- vortex in the electrode- 1, ob- 
tained by the similar maximization of equation (9) with 



Discussion of other unsustainable explanations 

Phase shifts in Josephson junctions may be induced by 
magnetic fields and currents in the junction area. As we 
discussed in the manuscript, we can rule out a possibility 
of direct penetration of the full vortex flux <& into the 
junction in a form of a Josephson fluxon[l]. Below alter- 
native plausible contributions are discussed in some more 
detail to show why they cannot explain the experimental 
observations either. 

First we emphasize that the measured flux-offset of 
I c (&) patterns is opposite to the trapped flux in the junc- 
tion. Indeed, if a positive flux (& v is induced by the vor- 
tex, then it will be necessary to apply the corresponding 
negative field to compensate this flux. Therefore the cen- 
tral maximum in I c (&) pattern, which correspond to zero 
total flux through the junction, would occur at a negative 
field, i.e, the whole I c (&) would be shifted by 



(10) 



Direct field of a vortex: Distribution of the magnetic in- 
duction in a bulk superconductor at distance r>A from 
the Abrikosov vortex is given by: 



B(r) 



$o exp(-r/A) 
2ttA 2 ,A7I : 



(11) 



If this field would stretch into the junction area, it would 
generate some Josephson phase shift, according to equa- 
tion (5). However, this field must have the same sign as 
the applied field, while in experiment we always observe 
that the effective trapped flux is opposite to the applied 
field. Thus, the phenomenon can not be due to direct 
field from the vortex. 

Effect of magnetic barrier: Barrier layers in our junc- 
tions contain Ni and are magnetic at low T in some cases. 
In the Nb/Pti-^Nlc/Nb junctions, the Ni concentration 
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was varied from about 67% (magnetic) , to pure Pt (non- 
magnetic). The magnetic barrier was used to reduce the 
critical current to a comfortable range below about 1 mA. 
The observed offset of the Fraunhofcr pattern is certainly 
not related to the magnetism of the barrier because a 
fcrromagnct would be magnetized in the direction of the 
applied field, in contrast with the experimental observa- 
tions. There is no reason for the magnetic moment of 
a ferromagnet to induce a $o/2 flux contribution either. 
Furthermore, no differences in behavior were observed 
for non-magnetic Nb/Pt/Nb junctions, also displaying 
flux shifts. The possibility to manipulate the state by 
transport current clearly demonstrates that vortices in 
electrodes are responsible for the observed flux offsets. 

Vortex currents: Distribution of the circulating current 
density around the Abrikosov vortex in a bulk supercon- 
ductor at distance r ^> A from the vortex is given by: 



Jq 



c$ exp(— r/A) 
'8^ 2 A3 ,A7a 



[1 + A/2r] 



(12) 



The influence of the vortex current on the Josephson 
phase shift in this case was considered theoretically in 
Ref. [4] . Here we provide similar estimations for the sim- 
plified case of a junction with a cylindrical electrode of a 
radius R with the Abrikosov vortex along the axis of the 
cylinder, z v = R, and the detector junction covering half 
of the surface of the cylinder. The Josephson phase shift 
is obtained by integration of equation (4) and the induced 
flux in the junction is obtained from equation (6): 

<f, if, Alfiv ~ s tZv exp (~ z v/ X ) 



2tt 



2A 2 V /^A 



-exp(-z t ,/A) 



\ 2zT a 



(13) 



The first positive term is due to the direct field of the 
vortex within the junction barrier, equation (11), and 
the second negative term is due to the vortex current 
at the junction interface, equation (12). If we substi- 
tute parameters for our mesoscopic Nb/PtNi/Nb junc- 
tions A ~ 150 nm for Nb, t = 20 nm and z v ~ 150 nm 
in equation (13), we obtain $> v /$o ~ 0.25 for the case 
of an anti-vortex. For the planar Nb/CuNi/Nb junctions 
the electrode thickness d — 70 nm is about two times 
smaller than A and the vortex is oriented perpendicu- 
lar to the thin superconducting film. In this case vortex 
currents decay slower (not exponentially, but quadrati- 
cally) with the distance from the vortex[3]. On the other 
hand the vortex is placed substantially further away from 
the junction z v > 300 nm. In any case we see that the 
contribution from the vortex current may be significant. 
However, it depends crucially on a number of parame- 
ters: z v , A, d, L. Since we observe the same behavior for 
two very different types of junctions with substantially 
different parameters, we conclude that the observed di- 
rect correspondence of the induced Josephson phase shift 



to the polar angle of the vortex can not be explained by 
the vortex current alone. 

Effect of stray fields: Stray fields from a vortex return- 
ing back through the junction area generate field in the 
direction opposite to the vortex field and would give rise 
to a flux shift in I c {<&) in the observed direction, simi- 
lar to the vortex current contribution. In fact, there is 
a well known fundamental connection between the mag- 
netic field outside the superconductor and the surface 
current in the superconductor. Therefore effects of the 
vortex stray field and the surface current at the junction 
interface are identical [follow for example the transition 
from equation (4) to equation (5)]. Discussion in terms of 
stray field lines from the vortex closing through the junc- 
tion area is perhaps more illustrative because it shows 
more clearly the nontrivial behavior of the trapped flux in 
the junction. Indeed, in the classical case, i.e., in the ab- 
sence of Josephson effect in the junction I c = 0, the stray 
field distribution should be governed by a simple magne- 
tostatics, which minimizes magnetic field energy outside 
the superconductor. In this case only a small portion 
of the total flux should go through the junction because 
field focusing in the narrow slit costs extra energy. To the 
contrary, we observe that a half of the flux A$ ~ $ /2 
is going through the junction for long enough junctions. 
Such quantization is not expected for the magnetostatic 
stray field. Most clearly the non-magnetostatic behavior 
of the trapped flux is demonstrated by the experiment 
with the stray field drain. As seen from Fig. lc, the 
drain opens an additional area for returning the stray 
field, which is more than twenty times larger than the 
area of the junction slit. It is clear that this should dra- 
matically reduce the magnetostatic flux in the junction 
slit. However, experimental data in Fig. 3 show that this 
is not the case. Apparently, magnetostatic stray field 
alone can not explain the observed phenomenon. 

Although the mechanism is still unclear, we believe 
that the presence of the Josephson junction itself plays 
a crucial role in the observed remarkable coincidence of 
the induced Josephson phase shift in the junction and 
the London gauge field of the vortex. Certainly, vortex 
stray fields and surface currents are involved. However, 
they are modified by the system to assume certain val- 
ues through the junction. In effect, the junction creates 
its own magnetic flux to complement any external stray 
fields to give a resulting phase shift of the observed mag- 
nitude. 
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